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1. INTRODUCTION 


The Agricultural and Resources Inventory Surveys Through Aerospace Remote 
Sensing (AgRISTARS) 1s a 6-year program of research, development, evaluation, 
and application of aerospace remote sensing for agricultural resources 
beginning In fiscal year (FY) 1980. The program is a cooperative effort of 
the National Aeronautics and Space Administration (NASA), the U.S. Departments 
of Agriculture, Commerce, and Interior (USDA, USDC, and USDI) and the 
U.S. Agency for International Development (AID). 

The goal of the program is to determine the usefulness, cost, and extent to 
which aerospace remote sensing data can be Integrated into existing or future 
USDA systems to improve the objectivity, reliability, timeliness, and adequacy 
of information required to carry out USDA missions. The overall approach 1s 
comprised of a balanced program of remote sensing research, development, and 
testing which addresses domestic resource management, as well as commodity 
production Information needs. 

The Technical Program is structured into eight major projects as follows: 

1. Early Warning/Crop Condition Assessment (EW/CCA) 

2. Foreign Commodity Production Forecasting (FCPF) 

3. Yield Model Development (YMD) 

4. Supporting Research (SR) 

5. Soil Moisture (SM) 

6. Domestic Crops and Land Cover (DCLC) 

7. Renewable Resources Inventory (RRI) 

8. Conservation and Pollution (C/P) 

The majority of these projects will make use of spectral information to 
estimate crop-related parameters. 


This literature review is Intended to be a comparative summary of the signi- 
ficant research results relative to the quantitative use of spectral data to 
measure the properties of leaves, plants, and canopies. The rapidly expanding 
interest in this topic has resulted in a spectacular proliferation of perti- 
nent literature. The authors have attempted to select a cross-section of the 
technical literature which illustrates the state of the art in remote sensing 
as applied to AgRISTARS research. This report is intended as a supplement to 
the recent commentary by Stuff and Barnett (1979). 

The dependent variables which have been the objective of most research and 
which will be treated in this paper have been divided into three categories: 

a. Yield and yield components [leaf area index (LAI), biomass, percentage of 
ground cover, and height] 

b. Stress factors (disease, insects, drought, malnutrition, physical damage, 
and salinity) 

c. Growth phenomena (emergence, leaf formation, reproductive phases, 
ripening, and harvest) 

Research has been conducted on several scales, including single leaves (and 
stacks thereof), plots, fields. Large Area Crop Inventory Experiment (LACIE) 
segments (5 by 6 nautical miles in size), and Landsat full frames (100 by 
100 nautical miles). Spectral properties have been measured by instruments 
which can be classified according to both wavelengths and number of wavebands 
used; e.g., continuous (narrow band), two-band, four-band, and thermal 
devices. Emphasis will be placed on the Landsat spectral range, 0.5 to 
1.2 micrometers (500 to 1200 nanometers). The following topics will be 
addressed: 

a. Fundamentals of Vegetative Reflectance - brief summaries of current 
opinions on value of measurements of the visable, near-infrared, and 
infrared water absorption bands 

b. Calibration: Corrections and Transforms - brief statements on principal 
procedures used to standardize spectral measurements over space and time 



c. Correlation Studies - discussion of representative or unique studies which 
have statistically related spectral measurements to yield and stress 
variables 

d. Temporal Trajectories and Crop Calendars - discussion of attempts to 
identify changes in plant appearance over time due to physiological 
development; i.e., growth phenomena 

e. Large Area Studies - qualitative evaluation of methods and performance of 
models or schemes for multiple site application 

f. Summary and Conclusions 

The following exclusionary rules have been applied in the selection of 

documents and titles to be included. 

a. Reject articles dealing primarily with classification. 

b. Omit articles and discussion of thermal and microwave techniques. 

c. Reject papers based purely on photographic interpretation. 

d. Omit discussion of radiative transfer models (atmospheric or plant). 

e. Reject articles not formally published. 



2. THE FUNDAMENTALS OF VEGETATIVE REFLECTANCE 


Eight comprehensive reviews on the fundamentals of vegetative reflectance have 
been identified. The authors are: 


a. 

Colwell, R. M., et al. (1963) 

e. 

Myers et al. (1970) 

b. 

Knipling (1967) 

f. 

Gausman (1974) 

c. 

Gates (1970) 

g* 

Myers et al. (1975) 

d. 

Knipling (1970) 

h. 

Harlan and Erickson 


A discussion of these reviews incorporated with views held 
this paper are presented in this section. 


(1977) 

by the authors of 


2.1 VISIBLE WAVELENGTHS (0.4 to 0.7 MICROMETERS) 

The visible region is the best understood part of the spectrum. This is to be 
expected since its interpretation can be tested directly against the millennia 
of agronomic experience. Spectrophotometric measurements of vegetation in 
this region are determined primarily by the amount of chlorophyll in the 
leaves. Normal concentrations cause most (80 to 90 percent) of the incident 
light to be absorbed. Because somewhat more absorption occurs in the blue 
spectrum (0.4 to 0.5 micrometers) and the red spectrum (0.6 to 0.7 micro- 
meters) than in the green spectrum (0.5 to 0.6 micrometers), the latter 
dominates vegetative reflectance patterns (signatures). As chlorophyll con- 
centration decreases from normal, the overall absorption in the visible region 
decreases but the green dominance remains. Hence, for individual leaves in 
the laboratory, an increase in albedo (overall visible reflectance) is usually 
the first visible indication of stress. However, in the field, many stress 
factors tend to reduce vertically observed leaf area (reduced canopy) leading 
to an increase in the ground contribution to spectral reflectance. Whether 
the overall visible albedo is increased or decreased depends upon the nature 
of the soil surface cover, including the degree of shadow and the moisture 
content. Row direction, of course, affects the degree of shadow. Also, the 
degree of contrast between soil and plant will affect the level of reflectance 
which is detected. This effect is especially important as resolution 
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decreases* Hence, Landsat albedo differences are particularly difficult to 
interpret. 

# 

Increasing stress normally leads to drastic reduction in chlorophyll content 
to the point where other plant pigments dominate spectral response. The 
carotenoids and anthocyanins thus become the main determinants if they are 
present, but some plants are low in such pigments and will simply increase in 
overall reflectance while band differences lessen. In other words, the spec- 
tral pattern becomes both lighter and brighter. Once again, the background 
soil surface contribution may have major effects, particularly at the Landsat 
scale of resolution. 

Finally, extreme stress leads from chlorosis (yellowness or blanching) to 
necrosis (dead tissue formation), producing a spectral signature that 
increases steadily with wavelength for individual leaves. This change is 
combined with an increasing and eventually dominant role of the background 
whatever it may be (green weeds, red soil, yellow debris). Clearly, the ner 
result is difficult to predict. 

In summary, stress effects in the visible range are well documented and under- 
stood, but the interpretation of an unknown target's signature can be exceed- 
ingly complex. The problem is further complicated by the fact that plants 
naturally display symptoms similar to those of increasing stress as they 
mature and ripen. Therefore, stress symptoms can be interpreted only if the 
growth history and stage of the crop are known. This point will be reiterated 
in the following discussion. One other complication that should be mentioned 
is that the dorsal (under) sides of leaves reflect more strongly in the 
visible range than do the ventral (upper) sides. Thus, wind-induced patterns 
can also confuse the interpretation of the spectral reflectance measurements. 

2.2 NEAR-INFRARED WAVELENGTHS (0.7 to 1.2 MICROMETERS) 

One of the most distinctive features of vegetation is that it reflects (and 
transmits) more in the near-infrared region than do most other natural 
objects. Many models have been proposed to explain the pathways of incoming 



radiation to and from vegetation which lead to the phenomenon. Most agree 
that this radiation Is first scattered and then reflected to the sensor by 
Index of refraction differences at cell walli and water Interfaces. 
Actually* the underlying cause of the strong reflectivity Is the almost total 
lack of absorption, which means that all Irradlance Is potentially available 
for reflectance or transmittance. This leads to considerable variability in 
reflectance, since It Is largely controlled by the varying leaf and canopy 
geometries of the multiple refractive Interfaces. (By contrast, the strong 
absorption by pigments In the visible band tends t/* stabilize reflectance pat- 
terns, so that only major physiological changes are Identifiable.) The sensi- 
tivity of the near-infrared means that this region Is a valuable source of 
Information when the variations In reflectance can be reliably associated with 
a particular physiological change. Much of the success of remote sensing In 
local areas has been In such cases. In fact, most of the specific examples 
discussed In the section of correlation studies involve the near- infrared. 
Unfortunately, Interpretation has proved to be difficult In large, complex 
areas where any of a number of factors may be affecting reflectance. Often, 
the directions of change of leaves and canopy have different signs. Also, 
canopy changes often lead to changes in the background exposure, which may 
further complicate the Interpretation. Several of the other caveats mentioned 
in the discussion of the visible range also apply to the near-infrared; e.g,, 
row direction and wind effects. 

2.3 NEAR-INFRARED l^AVELENGTHS (1.2 to 2.5 MICROMETERS) 

Laboratory studies Indicate that there Is very strong water absorption of 
irradlance at these wavelengths. However, the authors of this paper have not 
encountered many Instances in the literature where this fact has been tested 
for possible use in estimating water stress. Probably this Is due 1n part to 
the relative scarcity of Instruments which measure In these wavelengths. 

Since the Landsat thematic mapper includes provisions to collect data in these 
regions, their utility will probably be explored more thoroughly in the 
future. A few instances of such measurements are noted In the commentary on 
individual correlation studies. 


3. CALIBRATION! CORRECTIONS AND TRANSFORMS 


The Interpretation of Landsat measurements in terms of laboratory and field 
reflectance measure cents, the presumed research base, presents a serious 
problem because it is virtually impossible to establish a "maximum" from which 
to calculate true reflectance, In the absence of a calibration standard, 
Landsat "counts" have little meaning from day to day and from scene to scene. 
Possible solutions fall into three categories: 

a. Absolute calibration - This entails atmospheric radiative transfer models 
to correct for til significant distortion factors anc' thus permit 
estimates of reflected energy in absolute terms. 

b. Relative calibration - This requires an internal standard such as the 
scene mean, the scene maximum, and the scene minimum to permit conversion 
of "counts" to something like reflectance. 

c. Mathematical standardization - This relies on mathematical indexes, such 
as ratios and normalized differences to substitute for reflectance calcu- 
lations, There is no direct connection between these calculated values 
and reflectance. Hence, reflectance-based ground truth must be converted 
to these indexes, or a new ground-truth data base must be created using 
new measurements expressed in these mathematical terms. 

3.1 ABSOLUTE CALIBRATION 

The principal attempts at absolute calibration via atmospheric radiative 
transfer models are those of Lambeck (1977), Potter (1976), Potter and 
Mendlowitz (1975), and Turner et al. (1971). None of these has proved suffi- 
ciently useful to achieve general acceptance. However, the Lambeck model, 
XSTAR, is an integral part of the approach developed by the Environmental 
Research Institute of Michigan (ERIM). This approach will be treated 
separately. 



3.2 RgLATIVE CALIBRATION 

So far as the authors know, there have been no successful demonstrations of 
relative calibration methods although the scone minimum has been used In one 
version of the atmospheric correction model (ATCOR) developed by 
Potter (1976). 

3‘3 MATHEMATICAL STANDARDIZATION 

Most of the research aimed at normalization of Landsat data can be placed in 
the third category, that of mathematical standardization. The simplest of 
these transformations are the band ratios,. Of these, the most studied have 
been those of band 4 to band 5 and band 7 to band 5. The rationale for the 
first is that it is an indicator of the presence or absence of chlorophyll. 
The band 7/band 5 ratio is justified as a combined indicator of chlorophyll 
and mesophyll development. The reasoning is that it maximizes the vegetative 
nonvegetative difference by folding into one variable the facts that band S 
measures maximum chlorophyll absorption, whereas band 7 measures maximum leaf 
reflectance. (Some argue that maximum reflectance is band 6.) 

The use of band ratios implicitly corrects fcr major variations due to differ 
ences in Sun angle. For other transformations, a simple cosine correction is 
often used. In general, the digital data may be corrected to a reference 
sol or ele’'ation angle as follows. Let X, represent the Landsat signal in 
band i; the Sun-angle correction, X^-, is calculated as follnws: 

I cos On 

y a y. y 

''i cos 0 ^i 

where 0 is the solar zenith angle and 0q is the reference solar zenith angle. 
The resulting data will appear to have been acquired at the reference solar 
elevation angle. 





The second major type of transformation is the normalized difference, an 
example of which is the Transformed Vegetation Index (TVI) [Rouse 
et al., (1973)3. This transform has seen frequent use as an indicator of 
green biomass. The TVI may be calculated using either Landsat bands 7 and 5 
(TVI7) or bands 6 and 5 (TVI6). Let represent the Landsat signal in the 
band after Sun-angle correction, 


TVI6 « 


TVI7 = 


V 


Xg-X 



h “ ^5 


0.5 


0.5 


For some reason which we have been unable to identify, band 7 is not 
multiplied by two in this calculation, even though it would reduce the number 
of cases in which the 0.5 term is necessary to prevent taking th? square root 
of a negative number. Recent users [e.g.. Tucker, Elgin, and McMurtrey 
(1979); Tucker, Elgin, McMurtrey, and Fan (1979); and Tucker, Holben, et al . 
(1979)3 have dropped the use of the square root and are examining the 
normalized differences directly. 

The best known of the mathematical standardizations is probably the Kauth and 
Thomas (1976) Tassel led Cap transformation system used extensively by ERIM. 

The full flavor of this imaginative approach should be savored in the 
original. Briefly,- it consists of four mutually orthogonal transformations of 
the four Landsat bands, using empirically derived coefficients. Let X^ be the 
Landsat signal in the Hh band after Sun-angle correction and haze correction. 

Brightness = 0.433X4 + 0.632X5 + 0-586X6 + 0.642X7 

Greenness = -O.29OX4 - 0.562X5 + 0.600X5 + 0 . 491 X; 

The most important of the transformations is usually considered to be that 
which contrasts bands 4 and 5 with bands 6 and 7. Paradoxically, this trans- 
formation was christened Greenness even though it gives negative weights to 
bands 4 and 5, whose ratio forms the true basis of green. The other important 


transfonnation In this system is called Brightness, since It is essentially an 
albedo equivalent. The other two orthogonal transformations involve comple- 
mentary contrasts and are known as Yellowness and Nonesuch. They seem to have 
little informational content. However, Yellowness is incorporated in the 
XSTAR atmospheric haze correction algorithm. 


The XSTAR haze correction is applied after a Sun-angle correction. Let X* bo 

'A* 

a scene diagnostic signal for the ith Landsat band; let and X^. be coeffi- 
cients. Let Y represent the change in optical thickness from the reference 
condition. Let Y* be a reference yellow value for the scene. 


a s 

Y* B 

Calculate the following: 

a “ 

b = 

c = 

Y = 

where 



( 1.2680 \ /61.9\ 

1.0445 \ X* = I I 

0.9142 I I 83.2 1 

0.7734 / \33.9/ 

-11.2062 



-0.89952 X 
0.42830 I 
0.07592 I 
-0.04080 / 



The corrected va^ue* X^» is calculated from the Sun-angle-corrected band 
value, X^, as follows. 

X^ » e“''^(X^. - X.*) + X^* 

Misra and Wheeler (1977) and Wheeler et al. (1976) have derived a transform 
that is essentially equivalent to the Kauth-Thomas transform from theoretical 
considerations. 

The most common use of the Tassel led Cap transformation system is a pixel 
scatter diagram in Greenness (Y) and Brightness (X) space. The horizontal 
lower boundary of the plotted data is called the soil line. (Another defini- 
tion is the horizontal line below which falls 1 percent of the pixels.) Some 
investigators Ce.g., Thompson and Wehmanen (1977 and 1979)3 have subtracted 
the soil line value from Greenness and called the result the Green Number. 

The intent of this is to reduce the scene-to-scene variability in Greenness 
which is due to soil differences. Finally, an Idealized multispectral 
temporal Greenness and Brightness plot produces the pattern which suggested 
the name Tassel led Cap and is a form of spectral crop calendar. 

The Kauth-Thomas transform has been discussed at some length because it has a 
fundamental advantage over the other mathematical standardizations, such as 
ratios and normalized differences. This advantage is; the Kauth-Thomas 
transform preserves at least two dimensions in the data, whereas the others 
reduce the Landsat signals to a single index. We have already indicated that 
stress can only be measured if growth stage is known, and vice versa. 

Clearly, if only one independent predictor is available, the distinction 
between stress and growth stage effects will be difficult. There is an 
additional reason why at least two dimensions of variability are desirable. 

If vegetative cover is incomplete, the yield component percentage of ground 
cover is totally confounded with the other yield components, such as biomass 
and LAI. This is true because the percentage of ground cover is essentially a 
classification of vegetation versus bare soil. The other yield components are 
very highly correlated with the areal extent of the crop since bare ground 
possesses no plant attributes. Thus, plants can only be compared with each 


other after the soil contribution ■’’s removed; i.e., when complete coverage 
exists. Otherwise there is no way of distinguishing 4 metric tons of grass, 

50 centimeters high, and covering parts of the field equal to 1/2 hectare 
between the same amount of grass, 25 centimeters high, and spread evenly over 
an entire hectare. This confounding would cause serious problems for coeffi- 
cient estimation if a grass production model were attempted using spectrally 
determined percentage of ground cover 3 nd spectrally determined plant height 
as independent variables. The probl ;,n of confounding will always be a serious 
one. However, if it is ever to be overccme, the muUidimensionality of the 
data, whatever it may be, must be preserved. It should be obvious that it is 
undesirable to collapse all of the potential information into one variable and 
then expect to construct three models using that one variable to predict crop 
type and crop age and crop condition. 

In comparison, the ERIM approach is conceptually sound because it seeks to 
avoid these complications by first making an absolute calibration via XSTAR 
and then does a mathematical calibration which retains at least two dimensions 
of information. 


/3 


4. CORRELATION STUDIES 


The health status of a crop has always been the main concern of the producer 
because a crop stressed by nutrient deficiencies, drought, unfavorable 
weather, disease, insect infestation, and/or weed competition will net yield 
to its potential. Hence, the diagnosis of problems as early and quickly as 
possible allows managers to alleviate the stress and limit the yield 
reduction. 

Historically, this has been accomplished by field inspections and diagnosis of 
the symptoms, hopefully making correct judgments (e.g., distinguishing 
nitrogen deficiency from nematode damage) to permit proper action. This 
personal touch approach to crop protection requires a bit of luck and much 
time on the part of the inspector, who must be in the correct location at the 
proper time in order to spot problems while they are controllaole. 

With the advent of remote sensing, researchers recognized it as a useful tool 
for the rapid assessment of crop conditions over large areas. [See, for 
example, review articles by Heller (1978), Knipling (1967), and Wiegand, 
Gausman, and Allen (1972).] Even today people are impressed with aerial 
photographs and satellite images in which one can locate fields. Although 
some advancement in quantitative uses of these data have been made, the liter- 
ature continues to report studies in which stress is dealt with in qualitative 
terms. For example, healthy plants appear one color (or false-color) whereas 
stressed plants appear as another with no attempt to separate stresses caused 
by different agents or report the severity of damage. Few guidelines for 
application of their findings are offered [e.g., flyers et al. (1970), 
Malingreau (1978), Rohde (1971), Rohde and Olson (1971), and Seevers, Drew, 
and Carlson (1975)]. 

Fortunately for the producer, the number of quantitative studies that attempt 
the exploitation of information available from aerial photography and satel- 
lite imagery is increasing. To investigate relationships between spectral and 
agronomic properties (such as LAI, percentage of greesnness, percentage of 


ground cover, biomass, yield, disease, and nutritional status), researchers 
observe plants grovm either in normal production situations or under specific 
stresses, measure those properties of interest, and apply statistics to 
generate equations to predict agronomic variables from remotely sensed data. 
These studies will be discussed in two sections: crop condition and yield. 

4.1 CROP CONDIT ION 

As discussed earlier, the reflectance curves of normal and mal nurtured crops 
have a similar shape (e.g., where local maxima and minima occur), although the 
stressed crop may be more (or less) reflective in certain wavebands. 

Corn was raised under different nutrient treatments to determine whether 
spectral information should be used to differentiate between normal v nitrogen- 
deficient, potassium-deficient, and phosphorous-deficient plants (Younes et 
al., 1974). They reported that, dependent upon age, separation could be made 
between nitrogen-deficient, potassium-deficient, phosphorous-deficient, and 
normal plants in the visible bands. 

Al-Abbas et al . (1974) found potassium-, magnesium-, and calcium-deficient and 
normal corn to be more reflective than sulfur-deficient corn but less 
reflective than phosphorous-deficient corn between 0.75 and 2.6 micrometers. 
However, when they measured absorptivity, differences between treatments 
increased. 

According to Gausman, Escobar, and Rodriguez (1973), Mexican squash also 
reacted to nutrient stresses (nitrogen, phosphorous, potassium, iron, sulfur, 
magnesium, and calcium) by a shift to greater (or lesser) reflectance in 
certain portions of the spectra (0.500- to 0.900-micrometer range). Iron and 
magnesium deficiencies were distinguishable from those of potassium and sulfur 
in the 0.550- to 0.650-micrometer band. 

According to several authors, nitrogen deficiency increases leaf reflr'tance 
in the 0.500- to 0.700-micrometer band (as evidenced by the lighter to 
yellow appearance) in many crops: cotton (also as a result of sal in. 


reported by Thomas et al. (1966), cabbage reported by Thomas and Gerbermann 
(1977), and sweet pepper reported Thomas and Oerther (1972). In the latter, 
regression equations were developed at two wavelengths to estimate leaf 
nitrogen content from reflectance. 

Sorghum Is readily susceptible to Iron deficiency and becomes chlorotic as a 
result. Gausman, Cardenas, and Gerbermann (1973) took field and laboratory 
measurements of reflectance spectra and found that they differed from one 
another; although, In both sets, stressed plants were more reflective In the 
visible region than healthy plants. A later study by Gausman et al. (1975) of 
a sorghum field containing Iron-deficient areas confirmed the earlier findings 
of field reflectance. In addition, the 1975 study showed that chlorotic areas 
larger than 1.1332 hectares (2.8 acres) could be Identified from Earth 
Resources Technology Satellite (ERTS-1) data, using only one acquisition. 
Unfortunately, the investigators offered no discussion of areas which were 
identified as chlorotic but were actually healthy according to ground truth. 

Reflectance also changes with leaf water content because of various physio- 
logic responses to leaf water such as turgidity and transpiration (especially 
in the thermal bands). As with nutrient stresses, decreasing water content 
generally Increases the reflectance while retaining characteristic peaks and 
valleys in the curve. Thomas et al. (1971) found this to be the case in field 
and greenhouse studies of cotton, corn, and citrus. Linear regression allowed 
estimation of leaf water content from reflectance at two of the water 
absorption bands (1.45 micrometers and 1.93 micrometers). Cotton grown on 
saline soil was slightly more reflective (at 1.25, 1.45, and 1.93 micrometers)- 
than that grown on nonsaline soil. In addition, nonsaline-grown cotton was 
more turgid throughout the season. An earlier study also led by Thomas 
presented equations which related leaf moisture and reflectance In this region 
of the spectrum [1.45, 1.75, 1.93, and 2.22 micrometers (Thomas et al., 1966)]. 

The effects of the age of cotton and corn leaves on reflectance, water 
content, leaf structure, thickness, and chlorophyll content have been 
researched by the Soil and Water Conservation Research Division of the 


Agricultural Research Service in Weslaco, Texas. Gausman et al. (1971) noted 
that the age of cotton leaves caused the greatest variation in reflectance in 
the near-infrared (0.75 to 1.35 micrometers) with leaf reflectance increasing 
from node 2 to node 8, and then decreasing gradually downward to node 13. At 
all ages, maxima and minima in the reflectance spectra were retained. 

Gausman, Allen, Cardenas, and Richardson (1972 and 1973) confirmed the earlier 
cotton research and found similar effects in the reflectance spectra of corn 
of various ages. These experiments used field-grown crops. Not surprisingly, 
chlorophyll content at all ages was found to be significantly correlated with 
reflectance. Because of the continued compactness of corn mesophyll cells, 
reflectance in the 0.75- to 1.35-micrometer range changed little with age, 
whereas cotton leaves became more reflective in this region as the mesophyll 
cells expanded. Rate of nitrogen fertilization affects both chlorophyll and 
carotenoid concentrations, which in turn influence reflectance of visible 
light. Using five different nitrogen rates, Thomas and Gausman (1977) 
examined the relative effect of these pigments on the reflectance at 0.45, 
0.55, and 0.67 micrometers by cantaloupe, corn, cotton, cucumber, head 
lettuce, grain sorghum, spinach, and tobacco. They concluded that reflectance 
was primarily governed by chlorophyll concentration alone, and there was only 
a slight reduction in the total variance with the addition of carotenoid 
concentration. 

The state of the art in correlating remotely sensed data with diseased and 
insect-infested crops remains descriptive. However, the use of remotely 
sensed data in monitoring pests is feasible but requires a series of acquisi- 
tions to be accurate. Earlier workers sought previsual signs of stress, 
especially those due to insects and disease for the obvious reason of 
preventing outbreaks by controlling pests and disease in areas before they 
became unmanageable. The authors of this paper found no consistent evidence 
for previsual stress detection nor did Heller (1978) in his review of remote 
sensing applications. 

The Corn Blight Watch Experiment [MacDonald et al. (1972) and Kumar and Silva 
(1973)] was performed in the midwest during the 1971 growing season, the year 



following an epidemic of Southern Corn Leaf Blight In the Corn Belt. Biweekly 
aerial reconnaissance (Infrared and multi spectral) and ground truth were 
collected; fields were separated Into classes of blight severity. At early 
stages, blight was overestimated due to the (expected confusion with other 
diseases, Insects, drought, and malnutrition. As blight severity Increased, 
estimates from the aerial data became comparable to ground truth. 

Edwards et al. (1975) were able to distinguish between healthy orange trees 
and those stricken with young tree decline (YTD), using 16-band multi spectral 
scanner data from 457.2 meters (1500 feet). Tree health and YTD stage were 
assessed visually The best separation by computer Into healthy, slightly 
affected, or moderately to severely affected trees occurred In the 1.05- to 
1.1 -micrometer band. Classification was accomplished with an overall accuracy 
of 83 percent. When only two classes were considered (healthy versus all 
stages of YTD), classification attained 89-percent accuracy. (In this latter 
case, the best waveband for classification was 0.82 to 0.88 micrometers.) 


4.2 YIELD 

Most of the spectral yield models have been developed for large areas with 
Landsat data as Input and will be discussed in that section. 

One much published approach to yield estimates is the "stress degree day" 
models developed and reported in Idso, Jackson, and Reginato (1977); Idso, 
Reginato, and Jackson (1977); Idso, Hatfield et al . (1977); and Reginato, 

Idso, and Jackson (1978). The "stress degree day" 1s a summation of the crop- 
air temperature difference. Further inclusion of ground-based albedo measure- 
ments determines the minimum albedo which drives a linear regression which 
predicts the yield. The original model was developed for wheat in Arizona 
(both dryland and irrigated) with extensions to alfalfa and wheat in 
Cal ifornia. 

Tucker, Holben et al. (1979) modeled winter wheat yields from hand-.ield radi- 
ometer data (0.65 to 0.70 micrometers and 0.775 to 0.825 micrometers). They 
considered the spectral radiance of the two bands individually and in 


comblnatlcn (the ratios of infrared to red and the normalized lifference) over 
different periods (early-, mid-, and late-season) as predictors of final 
yield. Although the mid-season model gave high correlations, the regression^ 
relationships were not constant. Hence, Tucker et aU opted for integration 
of measurements over portions of the season as well as the entire season. 
Again, the midseason (stem elongation to anthesis) model had the greatest r^ 
(66 percent for the normal! zed-difference ratio). 

The reflectance of the short-grass prairie and its relationship to produc- 
tivity has been the continuing subject of research at Colorado State 
University and the USDA, directed by Miller and Pearson (1971), Tucker et al. 
(1973), Pearson and Miller (1973), Pearson et al. (1976), and Tucker in two 
reports (1977 and 1977). Several thousand measurements were taken to obtain 
relationships between reflectance and several agronomic variables including 
biomass, chlorophyll content, percentage of ground cover, and height. From 
these, they developed a hand-held radiometer to quickly and nondsstructively 
estimate biomass (from the ratio of reflectances at 0.80 micrometers to 0.68 
micrometers). Estimated biomass (as well as conventionally measured biomass) 
was used as ground truth for aircraft and satellite mapping of the condition 
of the short-grass prairie. Aircraft measurements taken were 12 band multi- 
spectral scanner, infrared photography, and Landsat thematic mapper. Areas 
were classified according to their aerially assessed conditions. Estimates of 
biomass were obtained from multiple linear regression of band values. 

In a report on grass canopy biomass classes, Tucker (1977) states that Landsat 
data could resolve only three classes of biomass on the short-grass prairie. 
Further, resolution was best at low biomass values (before full cover), 
implying that the classification was actually identifying only two levels of 
productivity. 

4-3 TEMPORAL TRAJECTORIES 

Temporal trajectories of spectral changes first proved their practical appli- 
cation in the classification of crop acreage. Since plants have basically 
similar reflectance curves, discrimination between crops is besV'done when the 





sequence (and timing) of colors Is followed as well as the reflectance at a 
particular Instant. 

Hlavka et al. (1978) reported a phenological crop classification scheme using 
mul t1 spectral , multitemporal Landsat data from five UCIB Intensive test sites 
in Kansas. They converted the data to a gray scale (0 to 31) and developed 
classification techniques by associating a growth state (determined in the 
ground truth) with each set of signals. CAlthough they did not define the 
growth states (up to 36 were used), they stated that, at any particular time, 
several states could be observed within a LAC IE segment and that results 
improved somewhat when states were pooled,] Growth-state signatures were not 
unique; hence, the earliest possible growth state was assigned to an obser- 
vation. These growth states had to occur In sequence for classification. The 
scheme Identified 89 percent of the corn and 83 percent of the wheat, with 
only about 5-percent false classification. They also reported that soil 
corrections and band ratio methods aided classification only slightly. 

In a crop, reflectance changes with time not only because of the changes at 
the cellular and leaf level but also because of the amount of soil cover. 

(Two common parameters of this are percentage of ground cover and LAI.) In 
addition, a series of acquisitions gives the most information on the progress 
and severity of stress on a crop. 

Other studies besides the previously discussed Soil and Water Conservation 
Research Division's cotton and corn projects and the Colorado State University 
and USDA's short-grass prairie projects have been conducted to determine the 
changes of reflectance with time and the reflectance relationship to biomass. 
Eleven crops (red fescue, quack grass, Kentucky blue grass, perennial rye 
grass, mixed grasses, potatoes, sugar beets, beans, barley, oats, and wheat) 
were grown In field plots and their reflectances measured several times 
throughout the season with a field spectrometer for purposes of classification 
(DeBoer et al., 1974). They found that the relationship between the biomass 
of grasses (in particular, Kentucky blue grass) and reflectance (at 0.871 and 
1.265 micrometers) was linear until the canopy closed; it became nonlinear 
after the canopy closed. 



Learner et al. (1978) observed seasonal reflectance patterns in two wheat 
cultivars in field plots with a ground-based speetroradiometer. They found 
that seeding rate affected the rate at which ground cover occurred but not the 
reflectance once 25-percent ground cover was reached. Wavelengths considered 
were midpoints of the four Landsat bands and in between the water absorption 
bands (1.65 and 2.2 micrometers). The latter were found to be good indicators 
of vegetative density, as soil is more reflective at these wavelengths and, 
consequently, easily separated from vegetative signals. Correlation coeffi- 
cients of reflectances with percentage of ground cover varied with the 
development of the crop, indicating the time dependency of reflectance. 

Kanemasu (1974) also considered field spectral measurements of wheat, sorghum, 
and soybeans on two soil types. He advocated the use of ratios as a means of 
removing solar el evation variations in order to compare data across the 
season. Although the ratio of reflectances at 0.54 micrometers to 0,65 micro- 
meters (corresponding to Landsat bands 4 and 5) could not discriminate between 
LAI and percentage of ground cover, equations were developed by Kanemasu to 
estimate these parameters from either the ratio or Infrared reflectance. 

Collins (1978) observed what he termed a "red-shift" of 0.007 to 0.010 micro- 
meter In the 0.745 to 0. 785-micrometer band In both sorghum and wheat at 
heading. In addition, the degree of heading and canopy density were 
Indicated, This red-shift was observable with high resolution spectral 
measurements taken of 300 fields during two airborne runs in California's 
Imperial Valley. A smaller red-shift was also noticed In alfalfa, cotton, and 
other nongraln’i. as they approached maturity. Collins hypothesized that the 
shift is due to increased chlorophyll concentrations during growth, causing 
stronger absorptance In the green band and the subsequent shift of the 
spectrum toward red. 

Tucker and others monitored the spectral changes of corn and soybean plots 
with time [Tucker, Elgin, and McMurtrey (1979); Tucker, Elgin, McMur^rey and 
Fan (1979); and Tucker (NASA in 79620, 1978)]. Red and Infrared reflectances 
and vegetative Indexes (from ground-based measurements) were correlated with 


agronomic variables (percentage of ground cover, plant height percentage of 
chlorosis, wet biomass, and dry biomass). All showed a definite time depend- 
ency. Tucker ei al. identified five distinct, spectrally measurable stages 
for both crops as follows; emergence to about 25-percent cover, rapid foliar 
growth and development, fall vegetative cover, onset of senescence and 
maturity, and crop maturity and harvest readiness. They comment on confusion 
due to stress and weeds and express hope that these considerations, along with 
stage information, can be incorporated into yield models. 

On a larger scale, the use of satellite data in the monitoring of the natural 
progression of vernal greening and autumnal browning across continents (termed 
"green wave" and "brown wave") was one of the first patterns to be noticed 
when satellite imagery became available. Consequently, the work done was more 
descriptive than practical, leaving countless regression equations which 
relate band ratios and vegetative indexes to the percentage of greenness of 
various target areas of known vegetation. Lilac, alfalfa, wheat, and range 
were considered in the cooperative study headed by Dethler (1974). Extensions 
of this work were published by his qoworkers [Ashley and Rea (1975) and Blair 
and Baumgardner (1977) on forests]. 

Carneggie et al. (1975) described the color change of California's annual 
grasslands with time as observed by Landsat. Ground spectral reflectance was 
measured at two sites to determine ratios to be used in interpreting the 
satellite data. They reported that the seasonal growth and dormancy cycles 
could readily be monitored. They found that the ratio of band 7 to band 5 was 
the best indicator of production and the green-feed period. These authors 
suggested that regressions or experience be useC to obtain quantity of produc- 
tion, acknowledging that, at the present, only relative range conditions could 
be deduced from Lf isat. 



5. LARGE AREA APPLICATIONS 


There are several problems in the evaluation of large-scale attempts to 
estimate stress from spectral measurements. In the first place, the ground 
truth compari son standard tends to be subject to considerable measurement 
error. Second, if bias is small, compensating prediction errors can indicate 
excellent agreement with ground truth overall, despite large standard errors. 
Third, some large area models have included nonspectral variables. This makes 
it difficult to evaluate the spectral contribution. Fourth, since failures 
are less apt to be published, the apparent significance of published results 
may be fortuitous. Some of all of these problems may be Involved In the 
consideration of the following examples. 

ERIM personnel [Colwell et al. (1977) and Nalepka et al. (1978)3 compared the 
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use of TVI, SQ75 (i.e., ^ mssb)» ^nd all-band models to predict winter wheat 
yields of several fields at three sites in Kansas. The study was based on 
previous results showing high correlation between percentage ’of ground cover 
and yield and between percentage of ground cover and SQ75. Results with all 
Landsat measures were remarkably good on both training and test fields. No 
clear superiority was demonstrated for any particular Landsat combination. 

The SQ75, at heading, was later used to predic\; whr-jt production in several 
counties (6 and 10) with excellent results. No formal evaluation of signifi- 
cance was performed. 

Schubert and Mack (1977) used percentage of ground cover to estimate yield. 

The authors call these biomass indexes. They actually are the percentage of 
pixels within an area of a particular vegetative class (wheat and all vegeta- 
tion) that are classified as having a closed green canopy based on the 
band 7/band 5 ratio. These indexes were calculated for three phenological 
phases: emerged (from emergence to 2 weeks before heading), heading (from 2 
weeks before to ?, weeks after heading), and ripe (wll subsequent dates up to 
harvest). When these indexes were compared with yields over 3 years (1973 
through 1975) and several sites, the correlation coefficients averaged above 



0.9 for the first two phases. The ripe phase correlations wer^ much poorer. 

A test of predicted versus published yields in a later year (1976) showed that 
only the wheat index was satisfactory. No formal statistical analysis was 
performed, but this type of variable appears promising, since it introduces 
spatial information. 

Helman et al. (1977) used a complex Landsat multichannel regression equation 
to estimate ' AI over a five-state winter wheat region. The LAI estimates were 
combined with an evapotranspiration model to estimate soil-water depletions 
and yields. The results were not evaluated statistically. 

Thompson and Wehmanen (1977, 1979) used a Green Index Number (GIN) technique 
developed in South Dakota to detect and monitor agricultural vegetative water 
stress over the U.S. Great Plains (USGP), the U.S.S.R., and Australia. The 
GIN is defined as the percent of pixels in a segment with green numbers 
greater than 15. Green number is obtained by subtracting the soil line from 
ERIM Greenness. The soil line is the lowest Greenness in the segment after 
eliminating the 228 (1 percent) pixels of lowest Greenness. A chi-square test 
showed that, in the USGP, the drought versus normal designations using GIN had 
a highly significant correlation with comparable splits based on the Crop 
Moisture Index (CMI) which is a measure of relative water availability. 

Results in the U.S.S.R. and Australia could not be evaluated statistically. 


6. CONCLUSIONS 


The technical literature demonstrates that while remarkable progress has been 
achieved in the uses of remote sensing in agricutural research a number of 
problem areas remain. The properties of electromagnetic waves are reasonably 
well known both in terms of propagation through realistic atmospheres and 
interaction with simple surfaces. However, the relationships between varia- 
tions in typical sensor data and events of agricultural interest are less 
clear. 

General applications of remote sensing to the estimation of plant growth, 
stress, and yield are dependent upon a number of factors. Improved resolution 
(spectral, temporal, and spatial) will certainly help. More immediate atten- 
tion can be given to a number of specific actions, which are feasible with 
existing sensors and within existing programs, these actions are: 

a. A strong and continuing commitment to a field measurements program is 
needed. This includes Landsat and aircraft sensor data coordinated with 
ground truth observations. Experimental sites should include normal 
agricultural as well as specially designed field experiments. 

b. Research data bases containing spectral, meteorological, and ground truth 
observations are a prime requirement. Wide temporal and geographical 
coverage for a number of crops is needed. Documentation and easy 
availability to the general research community are essential. 

c. Renewed research into signature extension techniques is needed for large 
area quantitative estimation of plant characteristics. The lack of con- 
clusive evidence that any current approach has general applicability to 
large areas argues that existing normalization techniques are inadequate. 
Priority should be given to the development of alternative approaches to 
this problem. 

d. Priority should be given to the development of spectral crop calendars. 
Yield components and stress indicators cannot be reliably estimated unless 
spectral crop calendars are available. 


e. Additional /•esearch is required on the effects of the averaging of 
spectral data at both pixel and field levels. The assumption that simple 
arithmetic averages are satisfactory is probably unwarranted. The use of 
spatial information (within field) should be investigated. 

f. Research into simple atmospheric correction algorithms should be 
supported. 
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